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Abstract 

We discuss precision measurements of the leading atmospheric parameters at a standard neu- 
trino factory with a detector that is sensitive to muons alone. The oscillation of the muon- and 
electron neutrinos in the neutrino factory beam to tau neutrinos adds to the muon events sample 
(both right sign and wrong sign) via leptonic decays of the taus produced through charge-current 
interactions in the detector. In particular, we study how this affects a precision measurement of 
the atmospheric mixing parameters and the deviation of <-> Vr mixing from maximality. In spite 
of the enhancement of the number of events due to the additional tau contribution, the determina- 
tion of the atmospheric mixing angle and the deviation from maximality will be poorer. We show 
that it is impossible to devise satisfactory cuts to remove this tau contamination. Neglect of these 
tau contributions will lead to an incorrect conclusion about the precision obtainable at a neutrino 
factory. 

PACS numbers: 14.60.Pq, 96.40.Tv, 95.55.Vj 
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I. INTRODUCTION 



Neutrino oscillations [1] are now well-established. The relevant mass-squared differences 
and mixing angles have also been studied in various experiments [2-6] ; the study of neutrino 
oscillation physics is now at the precision stage with many existing and proposed near-future 
experiments. A natural progression in the precision study is to probe oscillation physics with 
neutrino factory beams. 

In a neutrino factory, a /x"*" or /x~ beam from a storage ring decays down a long straight 
section of the beam-line that points at a far-end detector. The decay of the muons produces 
a spectrum of electron and muon neutrinos via /i+ Uf^e^Ue and ii~ — > i/^e~z7e- that is 
extremely well-understood. Charged- current (CC) interactions of the muon anti-neutrinos 
(neutrinos) in the detector lead to the production of /x"*" leptons. The electron neutrinos 
(anti- neutrinos), on the other hand, can oscillate into muon neutrinos (anti-neutrinos) which 
in CC interactions in the detector, produce ii~ (//"*") leptons, or muons with charge opposite 
to that from the unoscillated case and hence get detected as "wrong sign" (WS) muons. 

In general, the focus of studies [7, 8] with neutrino factories (and the many papers and 
effort towards specifying the parameters of such a future neutrino factory) has been to pin 
down the as-yet unknown neutrino oscillation parameters, in particular, the across-generation 
13 mixing angle ^13 and the CP violating phase Sep, in addition to other possibilities such 
as the mass hierarchy (or mass ordering in the 23 sector) and the octant of the 23 mixing 
angle 623 (whether it is larger or smaller than 7r/4). All these can be best studied through 
measuring the WS events, as discussed above, in detectors capable of charge identification. 
The so-called right-sign (RS) (having the same charge as unoscillated) muon events are useful 
for precision measurements of the leading atmospheric parameters 623 and the atmospheric 
mass squared difference, Am^, apart from their use in understanding cross-section and flux 
uncertainties. In particular, they are sensitive probes of whether i/^ ^ u^- mixing is maximal, 
when $23 — 7r/4; we shall refer to this value of 623 as being maximal. Such maximal mixing 
usually arises from an underlying symmetry [9] that would dictate the structure of the 
neutrino mixing matrix. Measurement of deviation from maximality is therefore of great 
signiflcance in developing models for neutrino masses and mixings. 

In this paper, we address the little-studied issue of contamination of the (right or wrong 
sign) muon events sample from oscillations of the muon or electron neutrinos or anti- neutrinos 
to tau neutrinos or anti-neutrinos, which, through CC interactions in the detector, result in 
tau leptons that subsequently decay to muons. While tau-production in neutrino-nucleus 
collisions has been studied in great detail, for example in Ref. [10], there are not many studies 
on the implications of the tau contamination (through leptonic decay) of the CC muon or 
electron events. 

Tau contributions naturally arise in multi-GeV CC neutrino-nucleus interactions. Al- 
though the tau CC events are kinematically suppressed due to the large tau mass, rrir = 1.78 
GeV, the total contribution from taus is substantial since <-> Ur oscillations are large, being 
driven by a nearly maximal 623- Here we focus on how this contribution affects a precision 
measurement of the deviation of ^ mixing from maximal, i.e., a precision measurement 
of 623 and its deviation from maximal. 

Tau contributions enhance both the RS and WS event rates in both the muon and electron 
sectors, since the leptonic tau decay fraction into both electrons and muons is about the 
same (^ 17%). In particular, we find that tau contributions substantially enhance the 
right- sign muon event rates, especially at small final lepton energies, independent of 613. 
However, due to a different parameter dependence, the tau contribution to muons worsens 
the determination of the deviation from maximality of 623 while leaving the determination of 
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|Am^| unchanged. We also show that it is practically impossible to devise satisfactory cuts 
to remove this tau contamination. Hence, studies (beyond those discussed here) at neutrino 
factories must first bring the relatively little-studied uncertainties due to tau production and 
decay under control in order to achieve the expected precision measurements with muon 
detection. 

The paper is organized as follows. In the next section, we briefly describe the physics 
of neutrino oscillations as is relevant for this study. We also list here the current limits 
on the oscillation parameters. In Section 3, we discuss the inputs such as the neutrino 
fluxes and cross-sections and preferred base-lengths. While we use standard packages for 
the deep inelastic scattering (DIS) cross-sections involving muon neutrinos, we use the mass- 
corrected leading order (LO) expressions for the tau contributions as discussed in Ref. [10], 
as well as comment on the next-to-leading order (NLO) contributions. In Section 4, we 
calculate numerically the event rates for various processes at an iron calorimeter detector 
capable of detecting muons with charge identification, such as the proposed ICAL/INO or 
the numerically mooted MIND. The numerical computation involves the calculation of the 
neutrino factory fluxes, the probability amplitudes using a solver for the neutrino evolution 
equation, which has been described in Ref. [11], calculation of the double differential cross 
section and in addition, the differential decay rate of the taus to muons, in case of the tau 
contribution. We conclude in Section 5. 



II. OSCILLATION PROBABILITIES 

With a long baseline, the large matter effects in the course of neutrino propagation help 
in improved sensitivity to many parameters. In general, the 3-fiavor probabilities, Pij, of 
flavor i/j oscillating into flavor i/j, depend on all the oscillation parameters (in vacuum): the 
mixing angles 6ij, i,j = 1, 2, 3, the CP violating phase, Scp and the mass-squared differences 
Aij = mf — mj, as well as the density of Earth's matter. The following analytic expres- 
sions [12] for these probabilities for neutrino propagation in earth matter in the constant 
density approximation, are useful to qualitatively exhibit the parameter dependence and for 
understanding the exact numerical results of the next section: 

Pe^ ^ sin2^23sin2 2^r3sin2A™ , 

^ cos^ 023 sin' 29^3 sin' A^^, , 

P;^ 1 - sin' 2^23 [sin' 6";;^ sin' A™ + cos' O'l^ sin' A™] - sin^ ^23 sin' 2^^ sin' A™ , 

^ sin' 2^23 [cos' sin' A^2 - cos' ^™ sin' sin' A^^ + sin' sin' A^J , (1) 

where terms involving solar mass terms have been ignored [13] and the superscript m refers 
to mixing angles and mass square differences in matter. 

For the atmospheric mass squared difference we use the notation [14] Am' = m' — |(m' + 
m|), so that a normal or inverted hierarchy is simply indicated by a sign (and not magnitude) 
change in this parameter. 

We use a Runge-Kutta solver to calculate the oscillation probabilities for various ener- 
gies and path lengths. Some technical details are given in Ref. [11]. All numerical results 
presented in this paper have been obtained using the density profile of the Earth as given 
by the Preliminary Reference Earth Model (prem) [15] and numerically evolving the fla- 
vor eigenstates through Earth's matter. In particular, the approximate expressions for the 
probabilities as shown in Eq. (1) are not used. 
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Parameter 


Best-fit value 


2-(T error 


A2i[10-^] eV2 


7 fi^+iJ-^^ 

r.DO_g 20 


7.25-8.11 


lAm^l [10-3] eV^ 


r, 40+0.12 


2.18-2.64 


sin2 012; [^12] 
sin^^ia; [013] 
sin^6'23; [6'23] 


0-304tn^6; [33.5°] 

n ni +0.016 

^•"-■■-o.oii 
0.501^6; [45.0°] 


0.27-0.35 
< 0.040; [11.5°] 
0.39-0.63 



TABLE I: Table showing currently accepted [16] best-fit values of oscillation parameters with 2a 
errors. In the case of the mixing angle ^13, a 2cr upper bound is shown. 

A. Existing constraints on the oscillation parameters 

The currently accepted [16] best-fit values of the oscillation parameters are summarized 
in Table I. The sign of the mass-squared difference Am^ as well as any possible deviation 
of ^23 from maximality (as well as its octant) are not yet determined. Also, there exists 
just an upper bound on the effective (13) mixing angle, while the CP violating phase bcp is 
unknown. 

Here, we will focus on precision measurement of deviation of ^23 from maximality, i.e., 
deviation from 7r/4. We base the detector near the magic base-line, at L = 7, 400 km, so that 
the results are insensitive to the unknown 5cp. We will assume that a combination of solar 

neutrino experiments and the KamLAND rcactor[2, 4] experiment completely determines 
the parameters A21 and ^^12. Since these play an insignificant role in determination of the 
atmospheric parameters, we set them to the best-fit values given in Table I and do not vary 
them in this analysis. 

B. Tau contributions and 023 

The z/^ ^ oscillations are large, being driven by a nearly maximal ^23. Furthermore, 
as seen from Eq. (1), P^.^ remains large even if 6*13 vanishes. Although there is kinematic 
suppression of the CC cross-section for tau production due to the large tau mass and the 
decay rate into electrons and muons is about 17%, there is still a sizeable tau production 
rate, which feeds into the muon production rate upon leptonic tau decay. In particular, this 
adds to the right sign muon event rate, and can even double it at low muon energies, since 
the tau decays preferentially produce low energy muons. It should therefore be taken into 
account in a precision analysis. 

On the other hand, as seen from the same equation. Per is proportional to powers of sin ^13 
and the tau contribution to WS muons is therefore small. It is of course comparable to the 
direct (from u^) wrong sign muon production rate to which it adds, although rather smaller 
since the large tau mass suppresses the cross-section; in addition, the contributing events are 
from further decay of the tau into muons and hence the rate is further suppressed. While 
the analysis of wrong sign muon events to determine parameters such as ^13, CP phase, etc., 
should perhaps include this contribution, it will have a negligible effect on the analysis at 
hand, viz., deviation of ^23 from maximality; however, in our numerical analysis we do include 
this contribution. 

In short, tau contribution to muon events from 1/^ oscillations is large and adds 

to the right sign events; its contribution from z/g ^ z/^ oscillations is small and adds to the 
wrong sign events. Hence there are additional muon events from Uj. in both RS and WS 
sectors to the direct events arising from (oscillated or unoscillated) z/^ interactions. 
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a. RS Events : The RS events are not very sensitive to the octant of 623. Notice 
that two of the terms of P^^ depend only on 2623 and are insensitive to the octant of ^^23- 
Further, for small ^^13, the octant dependent term in P^^ has negligible contribution, while 
P^T- is completely independent of the octant. Similarly, individually, both are sensitive to 
the deviation of $23 from maximality, but they have different dependences on this parameter. 

To leading order in 613, P^^ ~ 1 — sin^ 2623- F while P^t- ~ sin^ 2623- F, where F = sin^ A?g. 
Since the 623 dependent terms come with opposite sign, the combination of muons from direct 
production and from tau decays marginally decreases the sensitivity of the event rates to 
this angle (and its deviation from maximality). However, as we shall see later, it is almost 
impossible to remove the tau events. Any cuts that attempt to do so drastically reduce the 
direct muon events as well and hence worsen the sensitivity to the oscillation parameters. 
Thus, the combined contribution of events from direct muon production and those from tau 
decay is less sensitive to the deviation from maximality. Neglect of the tau contribution will 
lead to an incorrect conclusion about the precision possible for the deviation from maximality. 
This is the crux of our result. 

b. WS Events : Notice that Pg^ and Pe,- are both sensitive to the octant of the 23 mixing 
angle, 623, while their sum is less sensitive (indeed, the sum would be wholly independent of 
$23, if one were to neglect the differences in the i/^-nucleon and i/T--nucleon cross-sections). 
Hence the tau contribution affects precision measurements in the WS sector as well. 

Here we have focused on events with muons in the final state. However, taus also equally 
contribute to events with electrons in the final state. These events are interesting for other 
precision measurements and have been discussed in the context of i^^ Vf. appearance or the 
"platinum channel" where charge identification of electrons is also contemplated. Here the 
situation with respect to the tau contribution is drastically different. The tau contribution 
to electron events from oscillations (with subsequent decay of tau into electrons) is 

large and adds to the wrong sign events. Hence signatures in the "platinum channel" may 
suffer from these large additional contributions, independent of ^13. This platinum channel is 
well-motivated when ^13 is large (when the "direct" , that is, not tau-generated, contribution 
is large); in such a case, the tau contamination from z/g v-y oscillations also turns on and 
becomes a small addition to the right sign electron events. 

Hence the inclusion of tau leptonic decay boosts the signal for muon RS events and spoils 
the purity of the electron WS events. Details of this contribution to the electron events will 
be presented separately. 



III. THE INPUTS: FLUXES, KINEMATICS, CROSS-SECTIONS 

A. The neutrino factory fluxes 

We assume a basic muon storage ring configuration [17] with muon beam energy E\) = 25 
GeV, with 5 x 10^'^ useful decays per year. The muon antineutrino (neutrino) and electron 
neutrino (antineutrino) distributions for /j,^ decay in the muon rest frame are given by, 

d^ N- 2r^ 

oc =^[i3-2x)T{l-2x)P^cos9] , 



dxdfl in 
dxd^l An 



[{l-x)^{l-x)P^cose] , 



where P^ is the average muon polarization along the beam direction, E^ denotes the neutrino 
energy, x = 2Ey/mn and 9 is the angle between the neutrino momentum and muon spin 
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vector. 

For unpolarized muon beams, in the laboratory frame the boosted flux distributions are 
given by, 



Ely' (1 - /3 cos 99) 3m^ - AE^y (1 - /5 cos 





4n^ 


dydQiab 


ttL m 




24n^ 


dydOiab 


r 2 



^ E,y(i-/?cosv9) 



2^^1/(1-/? cos 99) 



where y — Ei,/Eh, 11^ is the number of useful muons per year, L is the base-line distance 
from source to detector, and 99 is the angle between the beam axis and the direction to the 
detector (which is assumed to be the forward or 2;-dircction). 

In general, the decaying muon beam has a finite divergence (~ mr or less) so that it has 
an angular spread about the central beam direction. Further, the neutrinos are emitted in 
a forward cone that becomes narrower as E^ increases. We can relate the neutrino-muon 
opening angle ip to the angles a (muon direction with respect to the z-axis) and {9' , (f)') 
(neutrino direction with respect to the ^-axis) as, 

cos </? = cos a cos 9' + sin a sin 9' cos 0' . 

The neutrino fluxes can be (analytically) integrated over the muon beam angle a, assuming 
a gaussian angular divergence of the muon beam around the z-axis with standard deviation 
taken to be [18] a = O.I/7, with the usual definition of 7 = Eb/nif^. The resulting neutrino 
flux is a function of {6', (f)') alone. A trivial integration over (j)' thus gives the neutrino fluxes 
as a function of the neutrino angle, 9' . The kinematic boost ensures that about half of the 
neutrinos are emitted within a cone 9' < l/'j. Even though the opening angle is very small, 
notice that the flux is not uniform in 9'; this can be important especially to understand the 
fluxes at near detectors as well as its spread at the detector. We average the forward flux over 
a small angle 9' < e = 0.3a where the intensity is fairly flat as a function of distance from 
the central axis to obtain the energy spectrum dN/ dE^ of the electron and muon neutrinos. 
This corresponds to knowing the neutrino beam direction to a precision of roughly 0.1 mr. 

The resulting neutrino spectrum with gaussian angular spread averaged over {9' < e, < 
(j)' < 2tt) is 

d^=£^7dQy^''\d^/^ = ^^f-'^^ -2(3-8,7)(l + cJe 

~y (7' - 1) [4 + Ce + + e-2-'3c, (1 + c,)] | , 



dE^ Eh Jdn 



\dydn/^ 



G 



-\y (7' - 1) [4 + c, + + e-2<^'3c, (1 + cj] I , (2) 



where c^, Sg refer to cos e, sin e. The (unoscillated) neutrino spectrum at a detector L = 7, 400 
km away is shown in Fig. 1. The neutrino and anti- neutrino spectrum of the same flavor 
(from muon beams of opposite charge) are the same if the muon mass is neglected, as has 
been assumed here. 

Note that, in the absence of oscillations, the number of events in the neutrino sector is 
always a factor of two or more larger than in the case of anti-neutrinos due to the larger 
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Ev(GeV) 

FIG. 1: Unoscillated electron and mu neutrino spectrum from a muon factory at a distance of 
L = 7, 400 km, with Ef, = 25 GeV muons. 

neutrino cross-section. Hence, studies of wrong-sign muons tend to focus on /i"*" beams where 
the wrong sign events correspond to fi~ from neutrino-nucleus interactions. Here, we are 
interested in the right sign events, and so would prefer fi~ factory beams; however, we sum 
over events from both fi~ and fi~^ beams, assuming an equal (and simultaneous) running 
time of 5 years in each case. 

B. The kinematics 

We are interested in processes that yield muons (of either charge) in the final state. 
These arise from CC interactions in which either muons or taus are produced, with the 
latter decaying into muons. Studies with neutrino factory beams typically use integrated 
cross-sections and aim at reconstructing the initial neutrino energy and direction; here we 
focus on the spectrum of the final state muon and hence require the detailed kinematics of 
the interactions, both for direct muon production and via tau decay. 

In the laboratory frame, a neutrino of flavor I = fi ot r and energy interacts with a 
nucleon of mass M and produces the corresponding charged lepton / at an angle 6i w.r.t. the 
incident neutrino direction, which we take to be the z-axis. In the case of Ur interactions, the 
azimuthal angle of the muon from tau decay, 0^, is also relevant: the tau is produced at a 
very forward angle while 0^ is restricted by the decay kinematics. The available phase space 
is restricted in both direct muon and tau production due to the constraint on the available 
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energy for the lepton: E_ < Ei < see Appendix A for the detailed expressions on the 
constraints. 

The effect of this pinching in available energy, for the case of a r lepton being produced, 
can be seen in Fig. 2 where the final hadronic mass mx is plotted as a function of Er- The 
notation is standard: m'j^- = W"^ = {p + q)"^, where p, q are the nucleon and the intermediate 
gauge boson 4-momenta in the laboratory frame. 




0.8 I ' ' ' ' ' ' ' 

23456789 10 

E,(GeV) 

FIG. 2: Kinematics of i/T-~nucleon CC interactions. The figure shows the allowed parabolas of 
constant cos Or in the mx-Er plane for = 10 GeV. The ends of the parabolas (at nix = 
nT'X,min = M) give the limits of the allowed tau lepton energy, which are severely restricted for 
cos Or away from 1. 

The allowed tau energy limits E± for a given cos 6r correspond to the end points of the 
parabola in the mx-Er plane. For example, for a typical neutrino energy E,, = 10 GeV, there 
is hardly any allowed energy region for cos 6*^ = 0.91. Hence, the tau leptons are produced in 
a very forward direction while the direct muons, due to their lighter mass, are less restricted. 

C. The cross-sections 

Since the energies of interest range from a few GeV to 25 GeV, the CC interactions include 
quasi-elastic (QE), resonance (Res) and deep inelastic (DIS) processes. 

Unlike many other calculations, we consider not the total but the double differential 
cross-sections for muon production with detailed energy and direction information so that 
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we have better control on the spreads in the muon energy and direction. Use of total cross- 
sections requires the reconstruction of the original neutrino kinematics and introduces larger 
errors due to large uncertainties in reconstruction of the associated hadron kinematics. The 
differential cross-section may be written in the form, 



where Gp is the Fermi constant, k = M^/i^Q"^ + M^) is the propagator factor with the W 
boson mass, Mw, Pi is the magnitude of three- momentum of the charged lepton and the Wi 
are structure functions corresponding to the general decomposition of the hadronic tensor 
for QE, Res and DIS processes. We have 



The detailed expressions for Wi are taken from Ref. [10] where the specific structure functions 
are hsted for QE, Res and DIS leading order (LO) processes. A few comments are in order. 

The Wi and W2 functions are those that appear in neutral current (NC) lepton-nucleon 
scattering processes as well, while W3 is a signature of the parity violating nature of weak 
interactions. The structure functions W4 and W5 only appear when the mass of the charged 
lepton is taken into account. For consistency, we use non-zero masses for both muon and 
tau leptons in the calculation; however, W4 and W5 are actually relevant only for the heavy 
tan. 

The quasi-elastic events dominate at lower energies. The separation between Res and DIS 
is arbitrary and determined by a cut, chosen to be > Went = 1-4 GeV. This still leaves 
= —q^ undetermined. However, reliable perturbative estimates for DIS can only be made 
when > few GeV^. A few percent of the events satisfy the cut on W for DIS but have 
very low (and hence very small values of the Bjorken scaling variable, x = Q'^/{2p ■ q)), 
thus making cross-section estimates for these cases extremely uncertain. Removing these 
events will lead to an underestimate of the actual number of DIS events; estimating them 
by scahng their to Q^j^ ~ 2 GeV^ overestimates the cross-sections since the structure 
functions are large when x is small. 

For both mu and tau interactions we use the CTEQ6 LO set of parton distribution 
functions. In this parameterization set, the starting value is Qq — 1.3 GeV, or the charm 
mass. Since backwards evolution is unstable, the problem lies with events with < Qq — 



The effect of either choice (removing the event or adding it by estimating its contribution 
at a larger Q^) is shown in Fig. 3 for both mu and tau interactions. In choice ci, for all 
events with < Qq, the parton distributions are estimated with Q^ = Ql, while as and x 
are kept at their true values. Hence this choice overestimates the event rate. In choice C3, 
all such events are thrown away, so this choice underestimates the event rate. 

The "true" cross-section would probably lie between the two extreme estimates and is 
probably close to the cross-section with choice C2, where events with Q'^ < 1 GeV^ are thrown 
away, while for events with 1 < Q^ < Qq, the parton distributions are computed with 
rescaled to Qq- 

Notice that the uncertainties are larger for the mu case (and also will be large for the 
electron case) and is less severe for the tau. 




(3) 





1.69 GeV^. 
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It should also be noted that this issue of perturbative tractability cannot be cured by going 
to higher orders in the calculation. For instance, the effect of including NLO contributions 
to the DIS cross-section in the Ci case in the limit of massless quarks and target hadrons 
is also shown in Fig. 3. The parton distribution functions are again from the CTEQ6 NLO 
set while the coefficient functions are given in Ref. [19]. While this is correct for the muon 
case, there are further mass corrections for tau interactions [20, 21] which are severe since 
they lead to effective higher twist terms. A discussion of these corrections can be found 
in Ref. [21]; however, the uncertainty from the small-Q^ events overwhelms the changes in 
going from LO to NLO. For the numerical calculations that follow, therefore, we use the LO 
expressions throughout, which are correct (including mass corrections) for both muon and 
tau interactions. 

Actual measurements at low and medium are essential to close this gap in our under- 
standing of neutrino-nucleon cross-sections. This may be possible in the near future with 
measurements from MINOS, Minerva, T2K, Argoneut, etc. [22]. 

In summary, we use the CTEQ6 LO parton distribution set throughout this paper, with 
choice Ci. The quasi-elastic, resonance and DIS total cross-sections for and Vj. (and 
their anti-particle)-nucleon CC interactions are shown in Fig. 4. The lower (upper) curve 
corresponds to interaction of the (anti)-neutrino with an isoscalar nucleon (average of proton 
and neutron). The actual numerical calculations use the differential cross-sections shown in 
Eq. (3). 




FIG. 3: Total deep inelastic scattering (DIS) CC i/T--nucleon cross-section as a function of the 
neutrino energy. The left figure shows the impact of different cuts on the cross-section. See text 
for details. The right figure shows the effect of including NLO effects to the c\ curve. The upper 
set of curves is for mu (or electron) DIS; the lower ones for tau. 
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FIG. 4: Total quasi-elastic, resonance and deep inelastic scattering (DIS) CC (anti)neutrino- 
nucleon cross-section as a function of the neutrino energy for and interactions. 



D. Tau decay 



Here we focus on the leptonic decay mode of the tau, for example u^-n^jy^, where the 

final muon from tau decay adds to the direct muon signal from i/^-nucleon CC interactions. 
While the branching fraction into muons (or electrons) is about 17%, we are again interested 
in the detailed kinematics of the final state muon and hence use the differential decay rates; 
see Appendix A for details. 

The typical decay rate for a tau with energy Er = 10 GeV is shown as a function of the 
muon energy and direction (with respect to that of the parent tau) in Fig. 5. Notice that the 
muon is emitted mostly in the direction of the tau (which was itself forwardly peaked) while 
its energy is typically highly degraded compared to its parent. Hence, muons produced in 
interactions mostly contribute to the forward event rate at low energy. The fractional decay 
rate has been plotted so that the area under the curve in both cases is 17%. 



IV. EVENT RATES IN A FAR-DETECTOR 
A. Preliminaries 

We assume that the neutrinos interact with a 50 kton iron detector such as the proposed 
INO/ICAL or MIND detector. Both /i"*" and /i~ beams are considered, with equal exposure 
for each. There are two kinds of contributions in the signal: 

1. The "direct" muons arising from CC interactions of muon neutrinos (or anti- neutrinos) 
on the target. 

2. The "tau- induced" muons arising from CC interactions of tau neutrinos (or anti- 
neutrinos) on the target, with the subsequent decay of the produced taus into muons. 
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FIG. 5: Tau decay rate as a function of the final muon energy E^^ and direction (inset) for tau with 
energy = 10 GeV. Here 9f^r is the opening angle between the muon and its parent tau. 

While the decay fraction is small, the tau events are large, as mentioned earlier, and 
hence this contribution is significant. 

In addition, there are right sign (RS) and wrong sign (WS) events for each of these. Table II 
clarifies the nature of different sources of muons in the final state. The charge of the RS/WS 
lepton depends on the charge of the muon in the storage ring. However, we do not identify 
the charge of the final charged lepton: all muon events are simply added. This is because 
the RS events are sensitive to deviations from maximality of 623] the WS events are very 
small compared to the RS events and their inclusion may only marginally worsen the results; 
however, the advantage gained in being "charge-blind" is significant, especially in the errors. 



Beam in 

Storage ring 


Neutrinos 
in beam 


RS 

Direct 


muons 

Tau 


WS 

Direct 


muons 

Tau 




T^e 




■ Pfj,r -(Jr-Tr 


iVe • Pefi ■ O'fi 


Ne ■ Per • O^r • 






■ Pf,f, ■ 


Nf, ■ P^r ■O't-Tt 


Ne ■ Pe^i ■ Cr^ 


A'e • Per ' CTt 'Tr 



TABLE II: The terms contributing to the production of right sign (RS) and wrong-sign (WS) 
muons from Ve and neutrinos and anti-neutrinos with /i^ beams in the storage ring. Negatively 
(positively) charged muons constitute the RS (WS) signal with neutrinos from a /x~ storage ring 
and vice versa with /U"*" beams. Event rates are proportional to the product of the oscillation prob- 
abilities, Pij, the neutrino fluxes, Ne^^, the differential CC cross-sections, (T^,t) and the differential 
decay rate, (in the case of r-induced muon production), as is symbolically shown. 
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The muon event rate is calculated as a function of the observed muon energy (and not 
that of the original neutrino that participated in the interaction). We have, generically, the 
number of muon events in a detector at a distance L, 



JE*)^" dEy Jo Je- di?^dcos6'^ 

n\,{E) = K / dE, 'IJ:' ' ■Pir(E,,L) / dcos^, / dEr- ^' 



dE'rdcos 9. 



r,- d£^^d cos 9fj,d(p 

where the superscript b refers to the beam-type, fi~ or /x"*" circulating in the storage ring; 
i — iJ,,e{jl,e) denotes the neutrino type produced in a beam. Here K = NtUy, 

Nt is the number of target nucleons (assumed isoscalar) and Uy is the number of years of 
data accumulation. The first expression in Eq. (5) corresponds to the production of direct 
mTion (denoted by the subscript D) production while the second corresponds to r production 
with subsequent decay into muons (denoted by r), with E being the observed energy of the muons. 
^max jg ^j^g maximum neutrino energy in the factory flux; the limits of integration are defined in 
Appendix A. Note that the integration over the final muon parameters, in case of tau production 
and decay, is restricted by an angular constraint; see Appendix A for details. The charge-blind 
events are obtained by summing over the index i and hence the total events for each beam are then 
obtained by the sum, 

T^ltiE) = J^inlniE) + nliE)) . (6) 

i 

Due to finite detector resolution, the observed lepton energy E may not correspond to the true 
lepton energy. We assume that the muons are reconstructed with a gaussian resolution width of 7% 
(o" = 0.07Efj_) for energetic, GcV muons. (We focus on the charged muon distribution and do not 
try to reconstruct the energy and direction of the original neutrino.) This identifies R in Eq. (5) 
as the Gaussian energy resolution function of width a. 

We assume a 90% reconstruction efficiency of muons which is charge-independent, although we 
implicitly assume an iron calorimeter detector with charge identification capability as in the MINOS 
[23] and the proposed ICAl/ino [24] or simulated MIND detectors. We use the final muon (from 
direct muon production or via tau decay) event rates accumulated over five years, to study the 
sensitivity to the deviation of the mixing angle 623 from maximality. 



B. The event rates 

At the magic [25] base-line, L = 7, 400 km, which is chosen here, there is no sensitivity to the 
CP phase 5c p which we henceforth set (arbitrarily) to zero. Typical event rates for the leading 
atmospheric parameters, Am^ = 2.4 x 10~^ eV^, ^23 = 42° (sin^ 2^23 = 0.9891), the reactor angle 
^13 = 1° (sin^ = 0.0003; sin^ 2^13 = 0.0012) and the solar parameters set to the central values 
given in Table I, are shown as a function of the observed charged lepton energy in Fig. 6. The panels 
show the direct muon production and tau decay contributions to the right sign (RS) and wrong 
sign (WS) observed muon events from fj,~ and beams. Since the neutrino-nucleon cross-sections 
are larger than those for anti-neutrinos, the /i~ ip'^) beams have larger RS (WS) events. Hence, 
/x"*" beams are preferred for studying WS events (where charge identification is crucial) which are 
sensitive to non-zero ^13, 23 mass hierarchy, octant of ^23 and the CP phase S, away from the 
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magic baseline. A fi~ beam with larger RS events would be desirable for precision of Am and ^23- 
However, using the full sample without charge identification is preferable, as the detection efficiency 
of the charge blind sample would be higher. Since the WS events are a small fraction of the large 
RS sample, using the sum does not worsen the precision of the parameters. We therefore add the 
events (RS+WS) from both fi^ and fj,~^ beams; the total event rate is also shown in the figure. 
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FIG. 6: Muon event rates over 5 years as a function of the observed muon energy, in bins of 2 
GeV at a 50 kton muon detector distant L = 7, 400 km away from a muon factory neutrino source. 
Right sign (RS) and wrong sign (WS) events from fj," and /x^ beams are shown in the upper 
panels. Contributions from direct muon production (denoted by D) and that of muons from tau 
decay (labeled as r) are separately shown. The left lower panel shows the sum D + t. Notice that 
there is a substantial contribution from tau decay to RS events, as seen in the right lower panel; 
also, the WS events have a negligible contribution to the total event rate. Oscillation parameters 
are as shown, with Am? in eV^ and other oscillation parameters as given in Table I. 
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It can be seen that there is a substantial contribution to the RS events from tau decay into 
muons. Since the tau Hfe-time is very short, these events will add to the observed RS events from 
direct muon production in f^-nucleon interactions. Since the tau contribution alters the sensitivity 
to the oscillation parameters as discussed in Section II B and illustrated in Fig. 8, we next discuss 
whether these events can be removed through various cuts on the final charged lepton energy or 
production angle. 

C. Cuts on tau contribution 

As stated earlier, tau production in neutrino-nucleon interactions is extremely forward-peaked; 
see Fig. 2. Furthermore, the muons produced in leptonic tau decay are also forward peaked; see 
Fig. 5. Hence, one obvious way to remove/reduce the tau contribution is an angular cut. Also, as 
can be seen from Fig. 5, the tau decay rate is large for muons produced with low energies, resulting 
in a substantial tau contribution to the muon events, at small observed muon energies (see Fig. 6). 
Hence, a muon energy cut can also be contemplated. 

The effect of cuts on the event rates is seen in Fig. 7. The panels show the effect of an increas- 
ingly drastic angular cut on the final state muons. It can be seen that the only cut effective in 
removing/reducing the tau contribution is one (with > 25°) that removes the signal itself! Alter- 
nately a muon energy cut, of about E > 10-15 GeV, can substantially remove the tau contribution, 
still leaving sufficient direct muons. However, such a large energy cut will will worsen the precision 
to which the mixing parameters can be measured as sensitivity is higher in the lower energy bins 
where matter effects are larger. In short, it is not feasible to cut out the tau contribution and still 
make a precision measurement, in this case, of the deviation of 623 from maximality. 

D. Effect of the tau contribution 

As stated earlier, the tau contribution alters the dependence on the mixing parameters, and thus 
alters the precision to which we can determine them. Fig. 8 shows the dependence of the direct 
and tau-induced muon events on ^23- While the tau events have less sensitivity to this parameter, 
the rate increases while the direct event rate decreases as 623 increases towards maximal 623 = 7r/4. 
However, the inclusion of muons from tau events alters the uncertainties considerably. 

A look back at the expressions in Eq. (5) clearly shows a major difference between direct and 
total (including tau contribution) muon production that is most easily appreciated by explicitly 
writing out the dominant RS contribution in each case: 

dn'j,{RS) ~ P^^dN^da^ ; 

dn\{RS) ~ P^^dA^^da^dF^ . (7) 

A near detector (where the oscillation probabilities are near zero and the survival probabilities are 
nearly 1), sensitive to muons, precisely measures the combination (dA^dci^), that is, flux times 
cross-section of the muons. This is what appears in the (RS) event rate for direct muon production 
and is therefore very well constrained from measurements at the near detector. Indeed, uncertainties 
are reduced to a factor of 10~^, mainly arising from differences in the shape of the flux at the near- 
and far-detectors. 

However, what appears in the (RS) event rate expression via tau production and decay is the 
combination dA^^do",-. While the flux is the same as for direct muon production and in fact very 
well understood, the same is not true of the cross-section. While the muon cross-sections have 
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FIG. 7: Effects of angular cuts on the tau contribution to muon events at neutrino factories. For 
details, see the text. It may not be possible to devise effective cuts to remove the tau contribution 
which must therefore be included in any analysis. 

large uncertainties inherent in any perturbative QCD calculation at low/medium Q^, these are 
compounded for the case of the heavy tau where mass corrections are large. Furthermore, since 
these contributions result from oscillations, no near detector can help reduce the uncertainties more. 
Hence overall uncertainties are much larger for the tau contribution than for direct muons and this 
exacerbates the problem of precision measurements with tau contamination. 

We now proceed to show this effect through specific numerical calculations. We use a combined 
overall normalization error of 0.1% for direct muon events (separate fiux and cross-section nor- 
malization errors are not required as for WS studies since the RS events are by far the dominant 
contribution and errors are kept well in check by studies at the near detector as described above). 
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FIG. 8: Variation of direct and tau induced muon events with 623- While increase in 623 decreases 
the direct muon events, those coming from production and decay of taus tend to increase marginahy 
with larger ^23- The conflicting behavior results in less sensitivity in the total events. 

We use a modest normalization error of 2% for the total (direct +tau) events; since the same (muon 
neutrino) flux contributes, this is essentially the error on the ratio of the tau to muon cross-sections. 
We use typical input values of (Am^, ^23, ^13) to estimate how well the generated "data" can be 
fitted, and calculate the resulting precision on the parameters. Again, we keep the solar parameters 
fixed at their best-fit values and set 5cp to zero. The best fits (and regions of confidence levels in 
parameter space) are obtained by minimizing the chi-squared according to the method of "pulls" 
[26] : 
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(8) 



Here the normalization uncertainty for the event rate in a bin is parameterized as a linear function 
of the pull: 

^Hn(e)=7^^^n(l + AiV0 , (9) 

where AA^ is either 0.1 or 2% for direct and total respectively and ^ is the pull that accounts 
for the overall systematic error on the rate. While the expression can be minimized over a set of 
pulls for each systematic error, this is not needed since the rate is dominated by the RS events as 
explained earlier. Here the theoretical rates in a bin correspond to a central value of = with 
^ = ±1 representing la errors. We follow the prescription of Ref. [26], where the chi-squared is 
first minimized over ^ and then minimized with respect to oscillation parameters to get the best-fit 
values and regions. 

E. Deviation of ^23 from maximality 



We present results for a typical sample set of input parameters, (Am , ^23) ^13 = 2.4 x 10 eV , 
41.9°, 1°). After minimizing over the pull, we minimize over Am^ and ^23; keeping the 13 mixing 
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angle fixed. Fig. 9 shows the allowed region in parameter space at 99% CL (Ax^ = 9.21). 

The solid line correspond to considering direct muon events alone and the dashed line correspond 
to using the total events, including those from tau decay. 
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FIG. 9: Allowed parameter space in space at 99% CL from CC events directly producing 

muons (solid line) and with the inclusion of muons from tau decay as well (dashed line) for input 
parameter values of (Am^, 6*23, 6*13 = 2.4 x lO'^ eV^ 41.9°, 1°). 

Note that the selected signal (all muons in the final state, independent of charge) is dominated 
by RS muons and is therefore insensitive to the octant of ^23 as well as the sign of Am^ (or the 
mass hierarchy). These can be probed only through a study of WS muons, as has been discussed 
elsewhere. Returning to the issue of deviation of ^23 from maximality, it can be seen that the 99% 
CL contour is much more constrained with direct muons than for total muons, including those from 
tau decay. In particular, it is the Am^ values that are smaller than the input (true) value, that 
broaden the contour and limit the discrimination. Hence the effect of adding in the tau contribution 
is a worsening in the precision with which ^23 and its deviation from maximality can be measured: a 
spread of ~ 2° if tau events are removed and ~ 4.5° when they are included. This occurs because of 
the conflicting dependence on ^23 of the two contributions, as has been discussed earlier. However, 
it can be seen from Fig. 9 that the inclusion of the tau events does not affect the determination of 
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(the modulus of) Am^ which is tightly constrained to better than 1% in either case. 

The largest true value of ^23 that can be discriminated from maximal through a study of muon 
events in neutrino factories is shown in Fig. 10, as a function of Am^. Again, ^13 is kept fixed at 
013 = 1°. While there is a distinct but mild dependence on Atv?, it is seen that tau contamination 
worsens the ability to discriminate 623 from maximal, thus making this measurement as well as an 
octant measurement harder than originally expected. 
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FIG. 10: The figure shows the largest (smallest) true value of 623 in the first (second) octant that 
can be discriminated from a maximal value of 623 = '7r/4 as a function of Am^ when only Direct 
(D) and total {D + r) events are considered. Here ^13 is kept fixed: ^13 = 1°. 



F. Effect of 6*13 

Due to the extraordinarily large RS signal, it is anticipated that some of the early measurements 
at a neutrino factory will be those of precision determination of Am^ and 623 and the deviation of 
the latter from maximality. Certainly a strong case for conventional neutrino factories of the type 
considered here will be made only when ^13 is small, below the reach of near-future super-beams or 
reactor searches. Hence we expect ^13 to be relatively unknown at the time of such measurements, 
typically restricted to a value ^13 < 2° (sin^ 2^13 < 0.005). The uncertainty in ^13 will worsen the 
allowed parameter space from both direct and total muon events. In particular, the dependence of 
the event rates (both direct and tau-induced) is such that values of ^13 larger than the true value 
lower the near maximal 623, and hence worsen the power of discriminating against it. In other 
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words, the tighter the upper bound on ^13, better will be the discrimination of ^23 from maximality. 
This can be seen from Fig. 11 where the chi-squared is plotted (for a fixed normalization, = 0, 
for clarity) for two different Am^ values, for the same input set as before. 



X 




40 42 44 46 48 50 42 44 46 48 50 



FIG. 11: Effect of ^13 on for input parameter values of {Am?, 623, 613 = 2.4 x 10~^ eV^, 41.9°, 
1°). Increase in ^13 lowers the near maximality, and hence worsens the power of discriminating 
against it, as shown for two typical values of Am^, for the same input parameter values. 

The effect again leads to a widening of the contour in regions where An? is smaller than the 
true value, as can be seen from Fig. 12. Here, the normalization is kept fixed to ^ = for clarity in 
understanding the effect of varying ^13 and we have assumed a bound ^13 < 2°, with the contour 
corresponding to allowed parameter values at 99% CL. Fig. 12 also shows, for completeness, the 
false minima in the second octant and in both octants with the inverted hierarchy. Similar results 
are obtained if the true 823 is in the second octant and/or the hierarchy is inverted, although the 
precision is marginally worse when the true 623 is in the second octant. 

Alternately, we can allow ^13 to vary freely, but include a prior on this parameter, through an 
additional term in the chi-squared to include the expected bound on sin^ 2^13 [27] from upcoming 
experiments, 

fsin^ 2^13 -sin^ 20*™^)^ 

x'^x' + ^ —2 (10) 

(^13 

where (T13 is the error on sin^ 261^;^^ which we take to be ai3 = 0.005. Using this modified definition 
of the chi-squared, and the same procedure as before, the 99% CL contours in the parameter space 
now correspond to Ax^ = 11.36 for 3 parameters, Am^, 623 and ^13. In Fig. 13 we show the largest 
true value of 623 that can be discriminated from maximal when both Am^ and ^13 are varied, as a 
function of the true Am?, given the normalization uncertainties as specified earlier. The inclusion 
of tau events still worsens the measurement, although the situation is not as acute as when ^13 was 
kept fixed. Again, there is better discrimination for larger Am^. 

c. Some miscellaneous remarks : We note that the ability to discriminate from maximality 
is marginally asymmetric between the two octants when ^13 is different from zero, as is the case 
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in all the examples shown. The results are better if the true value lies in the first octant but 
not significantly so. Furthermore, the discrimination capability improves as the uncertainty CJ13 
decreases. Finally, there is a non-trivial dependence of the results on the muon beam energy, Ef,. 
While the event rate increases with Eb, the number of "useful" events does not, since the sensitivity 
to mixing parameters is best for intermediate neutrino energies, Ei, ~ 5-10 GeV. This arises from 
the energy dependence of the oscillation probabilities as well as the energy dependence of Earth 
matter effects. It appears that muon beams with energies 15-25 GeV are best suited for such 
measurements. 
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FIG. 12: Allowed parameter space in space at 99% CL from all CC events producing 

muons, including tau events, with a fixed normalization = 0). The dashed lines corresponds to 
fits with a fixed ^13 = 1° and the solid lines result from varying ^13 with < ^13 < 2°. The input 
parameter values are (Am^,^235^i3 = 2.4 x 10~^ eV^, 42.0°, 1°). For completeness, the contours 
in the second octant as well as with the inverted hierarchy are shown. 
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FIG. 13: The largest (smallest) true value of 623 in the first (second) octant that can be discrim- 
inated from a maximal value of 623 = 7r/4 as a function of the true Am^ when uncertainties in 
Am^, and the normalization are taken into account. All the analyses use an input value of 
^13 = 1°. The solid (dashed) lines are for direct muon events only (total = direct + tau). The 
lighter obtained with an uncertainty of (T13 = 0.005 on sin^ 2^13 will approach the darker curves 
corresponding to fixed 6*13 as the uncertainty on ^13 decreases. 



V. CONCLUSION 

We have re-examined the precision to which the leading atmospheric parameters can be mea- 
sured at a standard neutrino factory. We have addressed the little-studied issue of contamination of 
the (right or wrong sign) muon events sample due to oscillations of the muon or electron neutrinos 
(anti-neutrinos) to tau neutrinos (anti- neutrinos). These, through charge current interactions in 
the detector result in tau leptons, which can then decay to muons, adding to the right as well as 
wrong sign muon events obtained directly, without tau production. The tau contribution worsens 
the precision to which 623 can be measured while leaving the determination of the atmospheric 
mass squared difference |Am^| unchanged. Hence the tau contamination worsens the ability to dis- 
criminate against maximal <^ Ur mixing. Any cuts imposed to remove the tau events drastically 
reduce the events from direct muon production as well and are hence impracticable. Uncertainties 
from this tau background have therefore first to be brought under control before making precision 
parameter measurements at neutrino factories. 
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Appendix A 



We detail the kinematics of neutrino-nucleus CC interactions with muons in the final state 
(direct muon events) and that of CC interactions with tau in the final state with subsequent tau 
decay into muons in the lab frame (tau-induced muon events). The events are characterized by the 
final muon energy and direction, i.e., the entire muon spectrum. Hence the detailed kinematics of 

the interactions, both for direct muon production and via tau decay, are required. 

We define the four-momenta of the incoming neutrino (k), target nucleon {p), produced r lepton 
{k') and the muon (gi) coming from the decay of the tau in the laboratory frame as, 

k = {E^,0,0,E^) , 

p = (M,0,0,0) , 

k' = {Er, Pt sin 0r, 0, Pt cos 9r) , 

qi = {E^,Pij,sme,j,cos(j)^,p^sme^sm(p^,Pfj,cos6^) , (11) 

where M is the (isoscalar) nucleon mass. The kinematics is standard and straightforward for the 
case of Ufj^ interactions (where there is no tau production, all components of k' are zero), with the 
angles of the direct muon produced in the CC interaction taking values < cos 6^ < 1 and (pfj, = 0. 

The threshold for the interaction is given by E^'^ = mi + m^/2M, for I = h,t; hence it is 
significant, Ei, > 3.5 GeV, in the case of tau production. Furthermore, in the case of I'r interactions, 
the tau is produced at a very forward angle, Cmin < cos Or < 1, in the lab frame, where 



M ml 



Finally, the angle (j)^ of the muon produced during tau decay is restricted by the decay kinematics 
to obey the constraint: 

2{E^Er - Pf,PT cos 0r COS 0^) - (m^ + ml) 

cos 6a > 1—7 — : — — . (13) 

2^?^p,-sm6'^sm6'^ 

In addition, the available phase space in both direct muon production and tau-induced muon 
production, restricts the available lepton energy to < Ei < E^, where I = ij,,t. We have, 

E±{E,,cos9i) = ^(6± V62-4ac) , (14) 

with 

a = {E^ + Mf -Elcos^Oi , 
b = {E^ + M){2ME^ + mf) , 

c = mfEl cos^ 01 + {ME^ + mf/2f . (15) 

At the limits, mx = mx,min = M, the nucleon mass. The notation is standard: mx = = 
{p + Q')^5 where q = k' — k is the 4-momentum of the intermediate gauge boson in the lab frame; 
mx is constrained to lie in the range ^/s — mi > mx > M, with s = {k + p)^ . 

The expression for semi-leptonic decay, r /J-eUe, of a tau into a final state muon, with four- 
momenta k' and gi respectively (defined in Eq. (11)), in the lab frame of the original neutrino- 
nucleon interaction, is given by. 



dT{Er,Ei^,0r,0i^,(j)^,) _ 2 Pti 
'^f 



dE^dcos9^d(l)i^ ^ ASErTT^^ 



3{ml + ml){k' ■ qi) - 4:{k' ■ qif - 2mlml] . (16) 
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The total decay width of the tau is = 1/cr where cr = ET-/mT{cTo) where the rest-frame hfetime 
is ctq = 87.2 microns. 
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